	v[0] = 0.25;
	q[0] = 3.0  * (y[1] - y[0]) * v[0];
	for(i = 1; i < L; i++)
	{
		v[i] = 1.0 / (4.0 - v[i - 1]);
		q[i] = (3.0 *(y[i]-y[i-1]) - q[i - 1]) * v[i];
	}
